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PREFACE 


The NSEG program was originally constructed by Mr. L. H. Leet of the 
United States Air Force Aeronautical Systems Division, Wright-Patterson 
Air Force Base. The code was subsequently modified and extended by 
Aerophysics Research Corporation under contract F33615-73-C-3039. The 
current version of the NSEG (Version III) extends the program appli- 
cability to higher speed (hypersonic turbo-ramjet) aircraft. It also 
includes various improvements generated by Mr. David T. Johnson of the 
Air Force Flight Dynamics Laboratory. The authors wish to extend their 
thanks to Mr. Walter Vahl of NASA for his extensive assistance during 
formulation and checkout of the turbo-ramjet propulsion system model now 
available in NSEG. The analytic basis of the turbo-ramjet model is 
due to Mr. Vahl. 

Mr. D. S. Hague of Aerophysics Research Corporation served as project 
leader for the present study. Dr. H. L. Rozendaal provided specialist 
support in the fields of propulsion system analysis and computer sciences 
Mr. R. T. Jones, formerly of Aerophysics Research Corporation, has also 
made significant contributions to the NSEG code in studies preceding the 
present one. 

Additional details and copies of the program deck can be obtained from 
NASA Langley Research Center. 
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NSEG - A SEGMENTED MISSION ANALYSIS PROGRAM FOR 
LOW AND HIGH SPEED AIRCRAFT 


D. S. Hague and H. L. Rozendaal 

AEROPHISICS RESEARCH CORPORATION 

SUMMARY 


Program NSEG is a rapid mission analysis code based on the use of approximate 
flight path equations of motion. Equation form varies with the segment type, 
for example, accelerations, climbs, cruises, descents, and decelerations. 
Realistic and detailed vehicle characteristics are specified in tabular form 
and a variety of layered atmosphere options are available. The mission 
specification is open-ended in that the upper limit on the number of flight 
segments to be included in a mission profile (currently one hundred and 
forty-nine) can be increased by increasing the size of'a single common block 
(C0MTAB) above its current size of 3000 words. The code contains an English 
language oriented input procedure for describing the mission segment sequence 
to be employed. In addition to its mission performance calculation capa- 
bilities the code also contains extensive flight envelope performance mapping 
capabilities. For example, rate-of-climb, turn rates, and energy maneuver- 
ability parameter values may be mapped in the Mach-altitude plane. Where 
suitable graphics capabilities exist these maps may be drawn by machine in 
the form of contour plots. 

The code contains several approximate flight path optimization capabilities 
based on Rutowski energy-like criteria. These flight path optimization 
formulations permit inclusion of minimum time or fuel flight segments 
and maximum range segments during climb or descent segments. Approximate 
take off and landing analyses are also performed. At high speeds centri- 
fugal lift effects are accounted for. Extensive turbojet and ramjet engine 
scaling procedures are incorporated in the code. Take off and landing 
analyses are also available which employ a high lift aerodynamic analysis 
model based on the Air Force Flight Dynamics Laboratory DATCOM method. 
Alternatively, user supplied high lift aerodynamics can be employed. 

This report is Volume I of 3 volumes. Total program documentation consists 
of: 

Volume I. Theoretical Development 

Volume II. Program User's Manual 

Volume III. Test Problems 



1. INTRODUCTION 


The NSEG code was originally developed by the Air Force and 
then modified under Air Force Contract F33615-71-C-1480 and subsequently 
extended to higher speed flight under the present study. Major changes 
made to NSEG are a complete code reorganization, more general vehicle 
specification, addition of energy maneuverability concepts, addition of 
plotting capability, development of new data input procedures, introduction 
of turbojet and ramjet scaling procedures and the accounting for centri- 
fugal lift effect. 

NSEG provides a generalized mission performance analysis capability based on 
approximate equations of motion for the state components.' 

{X i > * {V. h, Y, W, R, t> (1) 

In all flight modes the equations of motion 

{X.} = {^(V, h, Y, W, R, t; a, B A , N)} (2) 

are of an approximate nature. For example, in climbs, y is neglected. 
Approximate equations of motion are available for 


1. Take-off 

2. Acceleration 

3. Climb 

4. Cruise and loiter- 

5 . Descent 

6. Deceleration 

7. Landing 


Any number of mission segments may be pieced together to form a complete 
mission. Segments may be flown in either forward or reverse direction in 
any sequence specified by the user. This feature allows a direct solution 
to many two- or N-point boundary condition problems. 

A typical complete mission profile is illustrated in Figure (1).. The program 
may also be used to generate performance contour plots of the type illustrated 
in Figure (2) .. NSEG contains a variety of operating modes to aid in 
mission analysis which include 

1. Point performance characteristic evaluation where given 
{X} * {X^} , the function 

♦ = ♦(Xj) (3) 


is evaluated. 



2. Vector performance evaluation where given 

{X} ~ {»••••} Xi, X j , } j “ 1, 2, «•*, Nj (4) 

the vector 

{+} - Up = X i# Xj, ...} j = 1, 2, ...,Nj (5) 

is evaluated and the maximum or minimum value of <j> in the region 


% < x 5 < X J H 

is found by interpolation. That is. 


( 6 ) 


^ j = f(*-*> X^, X j , • ••) (7) 

3. Map performance evaluation where given 

m ±j = X i# Xj, ...} i = 1 , 2 , Ni 

j = 1, 2, Nj (8) 

the performance array 

[* ±j ] = [*(•••• Xi, X., ...)] (9) 

is evaluated over a rectangular mesh of points in the (X-^, Xj) plane 
and the resulting contours obtained in the manner of Figure (2) . 


4. Mission segment performance where given a state {X} , an approximate 
state equation, and a segment termination criteria, the state 
transformation, T^., which transform state i into state j according 
to 3 

{X}. T . -MX}, (10) 

1 J 

is accomplished. 


5. Mission performance where given a sequence of mission segments, the 
Successive state transformations 


{X> 1 -► {X> 2 + . . . + { X N _ 1 } -> {X N } (11) 

are completed. 


The analytic basis of program NSEG is presented below. Section 2 describes 
the propulsion system characterizations available. Aerodynamics are described 
in Section 3. Sections 4 and 5 present the vehicle weight model and the 
planetary representation. Section 6 describes the mission segment transformation 
calculation and Section 7 the mission segment options which use the transforma- 
tions. Program mapping capabilities are described in Section 8. 

Constants utilized in the NSEG code employ English Units (lb., ft., sec., naut. 
mi., and degrees Rankine) exclusively. It is imperative therefore that input 
quantities employ English Units; constants utilized in the illustrations 
reproduced in the text reflect this convention. Provisions of NASA Policy 
Directive (NPD 220.4) have been waived for those portions of this report that 
pertain to the NSEG computer code. 



2. PROPULSION 


2.1 Simplified Propulsion Systems 


All propulsive representations compute the vehicle fuel flow rate given a 
flight condition and the required thrust. Vehicle required thrust is computed 
internally by NSEG on the basis of instantaneous flight conditions. 

The maximum thrust, T max , is given by 


where 

^max " ^maxl * T max2* or T max3 

(12) 


T maxj“W M > h > 

(13) 

A throttle 

parameter, N, is determined by 



N * T ./T j = 1, 2, 3 

reqd maxj J ’ * 

(14) 

where T . 

reqd 

is the required thrust. Fuel flow is given by 



Wj = k • Wj (N, M, h) 

(15) 

or 

W 2 * k • W 2 (N, M, h) 

(16) 

or 

• • 



W. * k • W S CM, h) 

(17) 


The parameter k is a scalar for adjusting fuel flow to meet various specifi- 
cation requirements. Vehicle thrust, T, is determined by 


reqd 


(18) 


Within the program the simplified propulsion system input is scaled by 
appropriate factors to produce the following modified thrust and fuel flow 
data 


(T)! = T i /a ; i = 1, 2, 3 

(W)! = W.j/ClOO x 0 1 * 5 ) ; i = 1, 2 

(W)j = V q1 ‘ 5 


(19) 

( 20 ) 
( 21 ) 


where a is the atmospheric density ratio and 0 is the atmospheric 
temperature ratio. 
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In the dry and wet options (i = 1 and 2, respectively), the power setting is 
specified as a percentage of maximum; hence, the factor of 100 in the scaling 
equations for fuel flow in the cases i = 1 and 2. Again, in the dry and wet 
options, three power setting options are possible. These correspond to 


Maximum available in system 1 

Required to equal drag 2 


Specified value of power setting.... 3 


The corresponding fuel flows are determined from the power setting 


PS = FN/(FN Maximum Available) (22) 


It should he noted that power setting data is specified in the mission 
segment options of Section 6 and hence do not form part of the basic 
propulsion system input. 


2.2 Turbojet, Ramjet and/or Combined Engines 

with Inlet Precompression Effects (Combined Engine Option) 

The turbojet, ramjet, and /or combined engine simulation option provided in NSEG 
by the ENGINS subroutine provides a more realistic engine simulation option for 
high speed flight than those described in Section 2.1 due to the following 
features : 


1. Inlet flow field precompression effects due to airstream flow deflec- 
tion are accounted for in thrust calculations. 

2. A more accurate description of engine performance is provided by 
tabular engine data . 

3. Turbojet only, ramjet only, or turbo jet/ ram jet combination engine 
options are provided; in-flight option selection as a function of 
Mach number can be employed. 

4. Engine-related lift and drag forces are computed and accounted for in 
the net thrust and the desired vehicle lift. 

5. High altitude and Mach number real gas atmospheric properties related 
to engine performance are described by tabular data. 

The following sections provide an analytical description of the turbojet, 
ramjet and/or combined engine propulsion option and user-related information. 

2.2.1 Analytical Description 

The basic assumption of two-dimensional inlet flow is employed in the turbojet, 
ramjet, and/or combined engine simulation subroutine. The geometry and various 
flow regions relevant to the simulation are shown in Figure (3). The computations 
performed are divided into seven main categories as follows. 



X. Free stream properties at Station 0 which are a function of 
aircraft Mach number and altitude. 

2. Flow field conditions at Station 1, aft of the wing shock. 

3. Inlet computations which yield the turbojet inlet recovery ratio 
and ram drag. 

4. Flow field conditions at Station 2, aft of the engine inlet wedge 
shock. 

5. Turbojet thrust, airflow, and fuel flow computations, a function 
of flow field conditions at Stations 1 and 3 (at the turbojet 
compressor face) and throttle setting. 

6. Ramjet thrust and fuel flow computations, a function of flow 
field conditions at Station 2 and M^. 

7. Spillage drag computations, if any, are performed either in the 
turbojet iteration loop for the turbojet alone engine option or 
in the ramjet iteration loop for the ramjet alone engine option, 
or the combined engine option. Figure 4. 

Details of these computations are described in the following text. 

2.2.2 Freestream Properties 


For freestream conditions where M §3.5 and altitude exceeds 19812 m (65000 ft), 
real gas freestream total pressure (Pt(p and total temperature (Tt(p values are 
determined by the two-dimensional table lookup subroutine DISCOT. Linear 

interpolation is utilized. Plots typical of these data appear in Figure (5). 
Freestream total pressure in p.s.f. is then given by 

P t * 2116. 22*P* (P*in atmospheres) f2 3) 

If the above Mach-altitude conditions are not satisfied, total temperature 
and total pressure values are computed by first determining freestream 
static conditions (p 0 ,T ) from the 1962 Standard Atmosphere Subroutine 
ATMS62. Static enthalpy (H) and pressure ratio (p rQ ) are then determined 
from Table II, given the static temperature. Total enthalpy (H t ) can next 
be computed by the expression 

Hf - H ♦ (M 0 * » 0 ) 2 /S0073.2 (24) 


Given total enthalpy, the freestream total temperature and total pressure 
ratio (P r ) can be determined from Table II. Freestream total pressure 
is then to given by 


(F W * p o 


(25) 
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2.2.3 Flowfield Conditions at Station 1 


Given the freestream flowfield parameters at Station 0, their counterparts at Station 1 

aft of the wing shock are determined by the two-dimensional isentropic shock 

relations (NASA Report 1135) . These relationships are mechanized in the 

subroutine OBSHOK which performs oblique shock computations for wedge angles 

below the critical value and normal shock computations at higher wedge angle 

values. 


2.2.4 Inlet and Ram Force Computations 


Given flowfield parameters at Station 1 and the engine streamtube area, the 
ram drag is computed from the equation 

TURBOJET RAM FORCE = p x * * (M^ * a^ 2 + A^, ^ - P Q ) (26) 

J J 

where is the area of the stream tube of the airflow utilized by the 

* j 

turbojet. Ramjet ram forces are accounted for a priori in the ramjet specific 
impulse data. The ram force computed above, parallel to the vehicle wing 
under surface, is subsequently resolved into components in the lift and drag 
directions. 

The engine capture ratio (A c /Apy LL ) is determined from Figure 6, Table III, 
given M-p The air flow rate captured by the inlet then becomes 

W A c = W A puLL * ( - A c /A FULL' ) = P 1 A INLET * M 1 * a i (A c /A FULL' ) 

Capture ratio less than 1.0 indicates operation at a below design Mach number. 

As seen in Figure (3) the result of below design operation is that the 
wedge shock fails to intersect the inlet cowl resulting in inlet air 
spillage. This spillage produces forces (inlet spillage lift and drag) 
which are computed using a combination of continuity relationships and geometry 
as shown in Figure (4) . 

2.2.5 Flowfield Conditions at Station 2 

Given the flowfield parameters at Station 1, conditions at Station 2 are 
computed using subroutine OBSHOK as in 2.2.3 above. That is the two-dimensional 
flow behind the precompression (wing) surface shock is turned again through the 
inlet wedge angle (<5j) . 

2.2.6 Turbojet Airflow, Fuel Flow and Thrust Computations 


For a base size turbojet, airflow and fuel flow requirements at full throttle 
are first determined via the one-dimensional table lookup routine FTLUP. The 
corresponding maximum thrust is determined via the two-dimensional table lookup 
routine DISCOT. The corrected airflow (Wj0/<5) requirement is determined from 
Table V, given the total temperature at Station 3 (T t ) . This value is then 
used to determine the turbojet airflow at any flight condition by the equation: 



where 


(w A ) 

corrected 


* 6//T * F t 
'j 


(28) 


(W^) = the base engine corrected airflow rate, 

corrected 


« 



/2116.22, ratio of compressor face total pressure to 
sea level static pressure.' 


(29) 


F Tj 


. • T /518.67, 
r 3 


ratio of compressor face total temperature 
to sea level static temperature. 


(30) 


■ the turbojet scaling factor obtained by specifying 
(1) desired net thrust at sea level standard conditions 
or (2) desired turbojet airflow at sea level standard 
conditions 


Note: It is assumed that T t ^ = = T t (Figure 3), 


The fuel flow rate to air- flow ratio is determined from Figure 7, Table VII, 
given Tt • The maximum turbojet thrust to air flow rate is found from 
Table VI 3 using the two-dimensional table lookup routine DISCOT, given T t 
and the log to the base ten of the ratio of total pressure at Station 3 3 

at the compressor face and a reference turbojet discharge region static pressure 
P N which may either be the underwing pressure or the freestream static pressure. 
The above pressure ratio is obtained by first determining P* /P^ from Table IV, 
given M^, and multipyling by P-t /P^. ° 

The fuel flow rate, maximum thrust and specific fuel consumption are then 
computed as follows: 


( "fW ■ w (V“ A ) 

(31) 

f max • (f g /w a ) * (w A ) 

(32) 

(SFC W ■ CrW Wmax 

(33) 


Where the maximum net thrust of the turbojet given by the ' 

expression 

*MAX ** ^ am Forc ® 
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Note that the specific ^fuel consumption given above has units of (JLb. Fuel/ 
Sec)/Lb. Thrust. The non-standard time unit, is employed to be compatible 
with subsequent weight computations in NSEG. The fuel/air ratio Wp/W. is 
obtained from Table VII, given T . 

% 

The program is constructed to allow thrust required in the flight* direction 
to be determined either by a turbojet throttle setting input or a thrust required 
input. If the throttle setting option is exercised thrust required is computed as 


F = T * F. 
R J 


N 


C.35] 


where Tj is the turbojet throttle setting input. If Fg is input greater 
than the net thrust which the turbojet can provide in the combined engine 
mode, the excess thrust will be provided by the ramjet, if possible. 

Prior to entering the turbojet iteration loop the maximum turbojet net 
thrust is computed 


(F M ) ■ p uay - ^ DRAG 
N MAX MAX 


(36) 


and miscellaneous iteration parameters are initialized. 

The purpose of the turbojet iteration loop is to provide the net thrust 
required, if possible, totally by throttling the turbojet. This procedure 
requires that spillage drag, a function of turbojet throttle setting, be 
accounted for in the net thrust, i.e., 

F Tj " F R ’ C f TOt) Tj * C0S T - RAM DRAG - SPILLAGE DRAG (37) 


In the above expression y is the angle between the wing under surface 
and the freestream direction, RAM DRAG is the component of RAM FORCE 
in the freestream direction and SPILLAGE DRAG includes drag components 
due to inlet spillage and excess captured air spillage, if any. 

The logic within this loop first tests the turbojet throttle setting. 

If this value is less than one, iteration commences. Otherwise, turbo- 
jet throttle setting and total thrust are adjusted to their maximum 
values and a check for the combined engine option is made. If the combined 
engine option is not requested, an optional message will indicate that 
the thrust requested is higher than that which can be provided. Within 
the iteration loop, throttled turbojet performance parameters are computed 
by first determining the throttled engine SFC ratio from Table VIII, 
given the throttled net thrust to maximum net thrust ratio. The throttled 
SFC, fuel flow rate, and airflow rate then become simply 
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SFC « (SFC)^ • SFC/ (SFC) 


MAX C38) 

W T J ' SFC * (W T j (39) 

(Wa) tj - (Wf) TJ /(W f /W A ) C40) 

Since the throttled not thrust to maximum net thrust ratio required to 
compute turbojet airflow rate is a function of turbojet ram force which 
itself is a function of turbojet airflow rate, an iteration loop 
encompassing the above computations is required to ensure that the 
turbojet airflow rate used to compute the net thrust ratio is equal to 
that used to compute the ram drag. 

If as computed above is less than the previously computed inlet captured 
airflow rate, the excess is spilled by the turbojet and results in either 
a spillage drag, a ramjet thrust, or both. If is greater than the inlet 
captured airflow rate, the excess is assumed to be available from suck-in 
doors. 


If, (1) the turbojet throttle setting equals one and the combined engine 
option is selected or (2) a desired turbojet throttle setting has been 
specified by the subroutine input data the net turbojet thrust available 
in the freestream direction is computed by the expression: 


f Tj * (Ftot) t 

J 


• cosy - RAM DRAG - INLET SPILLAGE DRAG 


(41) 


and logic flow exits to ramjet calculations. 

If the turbojet throttle setting is less than one, the spilled air lift and 
drag components due to spillage of excess captured air (engine spillage) 
are computed as shown in Figure 4. The available turbojet thrust 

is updated to account for engine spillage drag by the equation 


F__ - (F , ) * cosy - RAM DRAG - INLET SPILLAGE DRAG (42) 

TJ TOT xj _ ENGINE SPILLAGE DRAG 


A maximum throttle setting check followed by a satisfactory thrust level 
check are performed; satisfaction of either results in a return to the 
calling routine. The iteration loop following no exit is a simple 
application of Newton's method which attempts to satisfy the required 
thrust level within one tenth of one per cent. A maximum of ten iterations 
is allowed to complete this requirement. The independent variable in 
this process is the total turbojet thrust. In the event convergence^ is not 
achieved in ten iterations, an error message is printed and a return to the 
calling program is executed with the engine parameters computed on the last 
iteration. 



2.2.7 Ramjet Airflow, Fuel Flow, and Thrust Computations 


The ramjet computation sequence is constructed similar to that for the 
turbojet. The maximum ramjet air flow rate to the full capture inlet 
air flow rate is computed from Table X given the total temperature at 
Station 2, T^. . The maximum ramjet air flow is then, simply, 
l 2 


( W A)RJ max * 


Wmax * 
Wfull 


Wfull 


C43) 


The ramjet actual fuel/air ratio divided by the stoichiometric 
fuel/air ratio , and ramjet specific impulse are determined from 
Table IX as a function of freestream Mach number. Engine spillage drag, 
if any, ramjet and total fuel flow rates, air flow rates, and net available 
thrust levels are computed in the ramjet iteration loop via the equations: 


0» £ ) m • * 'Wstoichiometric Ci4 > 


Stoichiometric 

F B (Wf) RJ * I SP* C0S Y “ SPILLAGE DRAG (45) 

(Wf) TO T * ( W ^RJ + (Wf) TJ C46) 

(VtoT - 0»a)rj + Wu C47) 

F NET " F RJ * F TJ (48) 


Where the ramjet air flow rate above is either due to the air flow 
spilled from the turbojet or the maximum ramjet air flow rate, which 
ever is less. In the latter case, the excess air is spilled at pressure 
?2 and engine spillage forces are computed. The term SPILLAGE DRAG 
above consists of only engine spillage drag in the combined engine option 
but includes inlet spillage drag in the ramjet only option. The ramjet 
iteration loop again applies Newton's method to satisfy, if possible, 
the thrust level requested within .1 per cent using the ramjet air 
flow rate as the independent variable. If this condition cannot be 
satisfied, an appropriate message is printed and a return to the calling 
routine is executed with the engine parameters last computed. 
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3. VEHICLE AERODYNAMIC REPRESENTATION 

All aerodynamic representations compute the vehicle drag given a flight 
condition and lift coefficient. Vehicle lift coefficient required is 
determined internally by NSEG on the basis of instantaneous fright 
conditions . 


3.1 Clean Aircraft 

Either of two aerodynamic representations may be employed for the clean 
aircraft as described below. 

3.1.1 General Form 
The clean aircraft drag is computed in the form 

C D " C Do * C Di t493 

where Cd o is the zero lift drag, and Cp^ is the induced drag. 

Both Cp o qnd Cp^ may be computed by three component summation; that is, 

CD 0 = Cp 0l + Cp 02 + Cp 03 (5.0) 

CDj = Cpij + CDi 2 + Cp i3 (51.) 

In Equation (50) each of the three component zero lift drags must be 
of the form 


Cp 0j = Cp 0j (h, M); j =1, 2, 3 (52) 

Similarly, in Equation (51) each induced drag component must be of the 

form 


c Di j * C Dij (C L , M); j - 1, 2, 3 


(53) 


3.1.2 Polynomial Form 

In this aerodynamic option the drag is computed in the component summation 
form 

Cp = C Dl + Cp 2 + Cp 3 (54) 

where 2 

c Dj - Cp 0 j + * 1 - • C L 2 + k2j(C L - Cl M in ) + ^3j C L 3 


j = 1, 2, 3 


(55) 
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and 


CD Oj 

= C Dq j (M) 

.(56) 

kl i 

■ Klj (M) 

(57) 

k2 i 

= k 2j (M) 

(58) 

l3 i 

= k3j(M) 

(59) 

:l min 

= Cl mtm (M) , the minimum drag lift 
M1N coefficient 

(60) 


3.2 Store and Pylon Drag 
Store and pylon drag is computed in the form 

Cog * Cpj + C^2 + Qlj (61) 

where 

CDj * Q)sj " N Sj + Q)gp_. * Nspj j = 1 > 2 , 3 (62) 

In Equation (7.2.25) the drag of a single type j store pair is 

c DSj = CD Sj (M) (63) 

The number of type j store pairs is Ngj . The drag of a single type j store 
pylon pair is 

c D SP j = c DsPj (M) (64) 

The number of type j store pylon pairs is NgPj 

3.3 Tank and Pylon Drag 
Tank and pylon drag is computed in the form 

C Dt - CDj + Cp 2 + Cd 3 (65) 

where 

c Dj = C D T .- * N Tj + c DTPj ' N TPj j = 2, 3 (66) 

In Equation (7.2.29) the drag of a single type j tank pair is 

CDtj * c DTj (67) 
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The number of type j tank pairs is Nj.. The drag of a single type j tank 
pylon pair is J 


CDTPj = ^Drpj (M) 

The number of type j tank pylon pairs is 


( 68 ) 


3.4 Centrifugal Lift 


Centrifugal lift accounts for the earth's curvature. In flight at constant 
altitude about a circular earth a vehicle is effectively falling to follow 
the earth's surface. In consequence the full vehicle weight is not supported 
during constant altitude cruise. It can be shown that the relieving 
centrifugal lift term is given by 


C 


L 

cent 


r 


W(V cos y)‘ 


R e +h 


qS 


G 0 Re 


4. VEHICLE MASS REPRESENTATION 


4.1 Overall Weight Empty 

The vehicle overall weight empty is given by the following equation: 
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W OWE = W BARE + 


2> 4 • "Si’fixed * 2>sPi ’ V’fi 
i=l i=l 


fixed 


where 

w 

OWE 

W BARE 

N Si 

WSi 

NSPi 

WSPi 

N? i 

w Ti 

Nip t 

W T p i 


' WT i ) fixed *S (NTP i ' fixed 

i=l 


( 69 ) 


i=l 

= overall weight empty 

= bare weight without stores, tanks, or pylons 
= number of store pairs type i 
= weight of one store pair type i 
= number of store pylon pairs type i 
= weight of one store pylon pair type i 
= number of tank pairs type i 
= weight of one tank pair type i 
= number of tank pylon pairs type i 
= weight of a tank pylon pair type i 


The suffix fixed indicates that only fixed tanks, stores or pylons which are 
not included in the payload must be included in the summations. 


4 . 2 Fuel Load 

The total vehicle fuel load is given by 

3 

W FUEL = W FINT + 2^ (N Ti * ^T^usable 
i=l 

where 

WpuEL = t 0 * 31 useable fuel weight 

W FINT = wei ght internal fuel 

= number of tank pairs type i 
w FTi = weight of fuel in one tank pair type i 
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The suffix usable indicates that the summation only extends over tank pairs 
which are not included in the payload. 


The total non-payload fuel on board at mission initiation is given by 


w FT n = Wp 


INT- 


i>Tt 


w , 


1=1 


FTi 3 * 5 


NPL 


( 71 ) 


where 

W FT 0 = initial fuel load 
Wp INT = initial internal fuel load 
= number of tank pairs type i 
Wpt^ = weight of fuel in one tank pair type i 

and the suffix NPL indicates the summation extends only over tanks which are 
not assigned to payload. It should be noted that the total fuel is specified 
directly by data input, and the internal fuel load is a computed quantity. 


4.3 Payload 

Hie total vehicle payload on board at mission initiation is given by 

3 3 

W PL ’ W PL int • "Sildrop SPi ' »SPi>drop 

i=l i=l 

3 3 

+ ^JN T . • W Ti ) drop + (NSTi * w TPi) drop C72) 

i=l i=l 

where 

Wpl = total payload 
WpL^t = total internal payload 

and the remaining quantities in Equation (72) are defined in Section 4.1. 


15 


' 5. PLANETARY REPRESENTATION 

A flat earth planetary model is employed. The gravitational force is 
a simple inverse square field. A layered atmosphere provides the following 
options: 

1. Tabular 1962 U. S. Standard Atmosphere 

2. Analytic 1962 U. S. Standard Atmosphere 

3. 1963 Patrick Air Force Base Atmosphere 

4. 1959 U. S. Standard Atmosphere 

5. January 1966 NASA Atmosphere, 30° North 

6. July 1966, NASA Atmosphere, 30° North 

7. Arbitrary Atmosphere Generated from Temperature and/or 
Pressure Variation with Altitude 

The N.SEG program basically computes a planar flight path. However, time 

to turn calculations are available; hence, a three-dimensional path can be 
analyzed by "folding" the path into a plane. Figure 8 . 


6 - FLIGHT PATH ANALYSIS 

Flight path analyses for take-off, climb, cruise, descent, and landing are 
included in NSEG. The analyses are all based on relatively rapid approx 
imate methods. Each flight path analysis model employed is described below 


6.1 Take-Off 


The take-off analysis performs the transfer 


{X} 


TO 


{X} 


50 


(73) 


where the suffix TO indicates state at beginning of take-off, and the suffix 
50 indicates state at the 15.24 m (50 ft.) obstacle. 


The take-off model was originally developed by Mr. Louis J. Williams of NASA's 
Advanced Concepts and Missions Division, OART, Ames Research Center. The 
program provides : 

1. Simplified high lift aerodynamics based on the USAF DATCOM 

2. A ground roll analysis 

3. Rotation logic 

4. Climb out to clear a 15.24 m (50 ft.) obstacle 

The take-off analysis is also available as a stand-alone code, program TOLAND. 
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6.1.1 Take-Off High Lift Aerodynamics 


The take-off model uses a self-contained aerodynamics package based 
primarily on the DATCOM methods. Angle of attack in the ground run and 
rotation maneuvers is determined from the vehicle geometry. In the ground 
roll 

(74) 


= a 


BG 


+ 


WB 


where 

etc = wing incidence in ground roll 
auG= body incidence in ground roll 
“WB® wing incidence relative to body 

In the rotated attitude 

“R ■ ““MAX - 1-0 * (75) 

The additional symbols are 
or = wing incidence following rotation 

a ®MAX = max; *- mum body rotation, usually determined by the tail dragging 
condition 


6. 1.1.1 Take-Off Lift and Drag 
6.1.1.1(a) Maximum Lift and Drag 
The wing maximum lift coefficient is given by 

Cl max = CCl max 5 base + ACl max + ACl flap (76 - ) 

with a corresponding angle of attack 

*’ ■ C “ Cl MAX> base (77) 

During take-off the maximum angle of attack, “max, is limited to 

“max ' o- 8 • *’ ( 7 «) 

In these two expressions 

Cl MAX =wing lift coefficient at the first peak, Figure 9 

^ L MAX^BASE =baSiC Wing maximum lift coefficient 

^Lmax “maximum lift coefficient increment due to taper and sweep 

A ^ l FLAP “maximum lift coefficient increment from flap deflection 

( aC L,MA^ BASE “basic wing angle of attack at maximum lift coefficient based on 
linear Cl q term. Nonlinear a increment is ignored. 
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The high lift aerodynamic model is a simplified DATCOM method for subsonic 
low aspect ratio, untwisted, symmetric section wings. Due to the low speeds 
encountered in take-off and landing, the DATCOM method is modified by the 
approximation 

8 =/l-M2 -1.0 (79) 

Clean wing contributions to Equations (76) and (78) are obtained from 
Figures 10 and 11. Figure 10 provides (Clmax^ BASE; Fig 111 © H gives AClmax* 
The wing taper ratio correction factors Cl and C2 of Figures 10 and 11 are 
obtained from Figure 12. In Figure 10, the take-off model is limited to the 
lowest curve, and the curve for M<0.2 is used in Figure 11. Angle of attack 
at maximum lift coefficient is obtained from Figure 13. 


Flap maximum lift coefficient increment is based on the expression 


Bf 


4 %lap - 10 - 95 < c i*>base 


0.5 


Bp 2 

4 1 


Cp 1/2 

(rM [.062S/F-. 000472 jF 2 ] 

C WE 


where 

^A^BASE 

BF 



Cf 


/ F 


= linear lift coefficient slope/degrees 
= flap span 

* exposed wing span 
= average flap chord 

= average exposed wing chord 

* flap deflection 


(80) 


6. 1.1.1(b) Ground Roll Lift and Drag 

During the ground roll, the lift coefficient is determined by 

CLg = 57.29 C Lo • sin(a G ) co»*(«e) * f cos (^j - ^ 

= F(a G ) sin 2 (a G ) cos(a G )/sin 2 (c« MAX ) (81) 

where CL a is the linear lift curve slope. The corresponding ground roll drag 
is taken as 

CD G - C Do + k • C Lg [C^ - Ci^ - AC LFLAp ] + C Dlq (82) 
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= zero lift drag coefficient 
= induced drag factor 

= lift coefficient at zero wing incidence 
= landing gear drag coefficient 

Rotation Lift and Drag 

The lift coefficient after rotation, Clr, is given by Equation (81) with a R 
replacing oq; that is, 

C Lr = p C«G) (83) 

The lift coefficient is subject to the condition that 

CLr 4 (1 • (84) 

This inequality constraint is imposed to prevent buffet or pitch-up problems. 

The drag coefficient after rotation is given by 

c Dr = C D 0 + kC LR[ C L R ~ c L q " AC L F LAp] + C Dlg (83) 

6.1.1.1(d) Lift and Drag at 15.24 m (50 feet) Obstacle 

The lift coefficient at a 15.24 m (50 feet) obstacle is based on the rotation lift 
coefficient. 

Cl 50 = c L r /( 1-!) 2 ( 86 ) 

The corresponding drag is given by 

C D S0 ■ C Do * k • Cl 50 [Cl 50 - C Lo - AC Lflap ] . C Dlg (87) 


where 

CD 0 

k 

CLo 

c dlg 


6.1.1-Cc) 


6.1.2 Ground Roll and Rotation 
The ground roll distance, Xg, is based on the expression 

13. 

X G ~ Ffo 

WTO " y 


■ _ f .G R-) 

G 2 l C L , J 


L R 


where 1 

C °GR = 2^ C °G + Cd R> 


( 88 ) 


(89) 
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and 

Ffo ■ take-off thrust 
Wto = take-off weight 

Uq' = vehicle rolling friction coefficient 
Time to reach the rotation point is given by 

T g = 1.1842 X G /V R 

where the velocity at rotation, Vr, is given by 


V R » 17 




( 90 ) 


(91) 


Rotation is assumed to occur instantaneously. 

6.1.3 Flight to Clear 15.24 m (50 feet) Obstacle 

The average drag coefficient between rotation and 15.24 m (50 feet) obstacle clearance 
points is assumed to be 


C °R50 = 2 (Cd R + C °50^ 

The distance covered in clearing the obstacle is given by 


Wtq 

50 - 0+2 - 745 %^ 


(92) 


(93) 


50 


H TO 


CD, 


R50, 


- 1 . 105 (777——) 


’TO C Lr 

Time to clear the obstacle after rotation is 

T 50 - V £1 -“M 11 V R> 


(94) 


Thus, total distance for take-off over 15.24 m (50 feet) obstacle is 

C G + X 50 


- X„ + X, 


The elapsed time is 


Total fuel used is 


T T0 * T G + T S0 


(95) 

(96) 
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Wp TO * F T0 ‘ T T0 /(IsP T0 ) 


(97) 



At the 15.24 m (50 feet) obstacle the flight path angle is obtained from 


where 


sin(Y 50 ) 



C D50 

CL50 


W 50 " W T0 ' Wp T0 

The corresponding rate of climb is given by 

RC S0 - 1.6889 V 50 sin(Y S0 ) 


(98) 

(99) 


( 100 ) 


6.2 Landing High Lift Aerodynamics 

The landing analysis closely follows the take-off analysis but in reverse 
sequence starting from the 15.24 m (50 feet) obstacle. The angle of attack at 
touchdown is 


TD BTD 
and in the subsequent ground roll 


1.0 + 


“WB 


where 


°LR = a BLR + °WB 


(Xpjj = wing incidence at touch down 
“btd® bod y incidence at touch down 
°LR “ incidence during landing ground roll 

body incidence in landing ground roll 


( 101 ) 


( 102 ) 


6.2.1 Landing Lift and Drag 

The wing maximum lift coefficient during landing and the corresponding angle 
of attack are given by Equations (76) to (78) . Flap incremental lift is 
given by Equation (80) . It should be noted that the landing configuration 
parameters such as flap angle and permissible body angle of attack will 
normally differ significantly between the take-off and landing configurations. 
At the 15.24 m (50 feet) obstacle, configuration lift is assumed to|be 


c M.so " tl0!) 

At touchdown, Cl^, is based on Equation (80) using an)> that is 

Cl TD * p (°TD) + AC Lpyu> (104) 
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where AClrlap is given by Equation (80) using the landing flap setting. 

The inequality 

c Ltd < C 105 ) 

is used. 

Similarly, during the subsequent landing ground roll, 

Cl L r = F ^ a LR^ C 106 ) 

Drag coefficient at the 15.24 m (50 feet) obstacle, Cdlso> is given by Equation (82) 
using appropriate landing coefficients. Drag at touchdown, Cdtd> given 

by Equation (7.1.8) using touchdown coefficients. Drag during the landing 
ground roll is given by 

C °LR “ C °o + k * C Llr^ C l LR " <V + C °LG + Cd CHUT (107) 


where 


CDCHUT = lan<i i n g parachute drag 
All other symbols are defined in Section 6.1.1. 


6,2.2 Flight from 15.24 in 
Velocity at touchdown is assumed to be 

V TD = 17,16 V S, C Ltd 
T he ground distance covered from 15.24 m 


(50 feet) Obstacle to Touchdown 

(108) 

(50 feet) obstacle to touchdown is 


50.0 + 2.745 ( g ^ ) 
m S>Cl td 

L50 ^TD50 Ft, 


(109) 


where 

Xlso * flight distance from 15.24 m (50 feet) obstacle to touchdown 

Wl = landing weight 

Cd TD 50 * 7( c Dso + c D to )» the average drag coefficient 

Fl * approach thrust 


Z2 



Rate of sink at the 15.24 m (50 feet) obstacle is 


__ , lr r , CP L50. ,*150, ,, 

Rs 50 * 1,69 ■ C W L 


Flight path angle at the i 5.24 m (50 feet) is given by 

r C °L50 % , F L50^ 


sin(Y L50 ) * C 


c Llso w l j 


The ground roll distance is given by 

„ »- ot <ctbs> 


GL 


U+2<CD LG -U • c Llr) / C LxD 


Total landing distance is 


X L * X L50 + X GL 


(110) 


( 111 ) 


CH2) 


(113) 


6.3 Alternative High Lift Aerodynamics 


Where detailed high lift aerodynamic values are known the DATCOM calcula- 
tions may be bypassed. In this case the values of Cl > Cl and a MAX 
(CLAT0T, CLMXT0, ALMXT0) are specified directly in “ MAX 
program ' input . 

6.3.1 Acceleration at Constant Altitude 


The level flight acceleration segment performs the operation 

; X2 

(X) 2 * {X> 1 + / {X}dM * {Xj} + £{AX). 

where {AX}^ is the state change in accelerating from M to M + 


(114) 


AM. 


Given {X}^ ■ {V, h, y, W, R, t}^, Tj and Dj, then the velocity change is 
V * a s AM ( 115 ) 
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where a. is the speed of sound at the acceleration altitude and 


*i ‘ *» iT&i 

The approximate time to accelerate from to Mj +AM is 


At* = AV/Vi 

The corresponding approximate weight change is 

AW* = W. At' 

1 


and 

and to the first order 


W i+ i = W. - AW 


Vi " «o fT 3 


C116) 

Cl 17) 

(118) 

(119) 

( 120 ) 


i+1 


W i+1 can be obtained at the new Mach number. The mean acceleration is now 


5 - 2< V i * W 

The revised time to accelerate is 
At = AV/V 

which gives a weight change of 

SW = }(»! * * ltl )M 
and a range increment 

4R » 1 tV £ . V ltl )lt 

The state incremental vector {AX}^ is therefore given by 


(121) 

( 122 ) 

(123) 

(124) 


~AV~ 

= 

a s AM 

Ah 


0 

AY 


0 

AR 


y(Vi + V i+1 ) At 

AW 


j(Wi + Wi+i) At, 

. At. 


2AV 


(Vv i+ i) 


(125) 
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6.4 


Accelerating Climbs 


All accelerating climb paths are formed by a sequence of elemental straight 
line arcs in the Mach-altitude plane. On any arc the vehicle flies from 
(M^, h^) to (Mi+j, hi+i) . Since the vehicle is climbing 


h. - > h. 
i+I 1 


(126) 


The typical arc. for a climb path is shown below.. The Mach-altitude plane 



can be transformed into the velocity-altitude plane as follows: 



V = V(h, 

M) 


(127) 

so that 

AV = iY. ; 
AV 3h 

6h + 

s- 

(128) 

or 

dV 3 V 

3V 

3M 



dh " 3h 

3M 

3h 

(129) 


^ + 
- 3h 

dM 
a dh 


(130) 

where a is 

the local speed of 

sound 



V = aM 



(131) 


Now 3V/3h is the change in velocity with altitude at constant Mach number, 
and from Equation (131) with M constant 


3V 

3h 





£Tr 

dh 


(132) 
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where T R is temperature ratio, T/Tsl> so that 

9V _ V 9a_ dTR 
9h V * 9T R ' dh 

and from the atmospheric model 

1/2 

a = 1116. 45 (T r ) 

9a a 
9T R = 2Tr 

Substituting into Equation (133) 

9V _V_ ' dT R 
9h 2 Tr dh 

Substituting Equation (136) into Equation (130) 


(133) 

(134) 

(135) 


(136) 


dV V dTR dM 
dh 2 T r dh a dh 


(137) 


Equation (137) is used to define the required variation of velocity with 
altitude over an elemental climbing arc. 


Now the rate of climb is 


dh 

dt 


RC 


(138) 


or 

dh - A+ 

RC “ dt (139) 

Assuming rate of climb varies linearly with altitude in the elemental arc 


RC = a + bh 

Substituting into Equation (139) and integrating 


(140) 


or 


hi 


/ 


h2 


dh 


(a+bh) 



h+i 


t. 

x 




h. ,-h. n 
i+l i 

RC i+l _RC i 



(141) 


(142) 
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where 


a = RC^ 


(143) 


b = RCj+i-RCj ( 1 44) 

h i+l~ h i 

The vehicle rate of climb is computed under the assumption that thrust is 
aligned along the velocity vector as shown below. 



mV = (T - D) - W sinY (145) 


but 


Combining 


mV 

Equations 


dV dh W dV 

= m dh ' dt g dh 

(145) and (146) 


V sinY 


(146) 


and 

so that 


siny 


T-D 


W r IdV 
" l g dh 


+ 1 . 0 ] 


cosY = /T-sinZV 



(147) 


(148) 


(149) 
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Equation (149) can be evaluated at each end of the elemental arc to obtain 
RCi+i and Rq. Hence, At, the time to traverse the elemental arc, is given 
by Equation (143) . Similarly, the flight path angle at each end of the 
arc can be obtained from Equation (147). It should be noted that if siny. 
Equation (148), is greater than 1.0, the approximate climb analysis is 
invalid. If this condition occurs, the thrust is reduced to produce a 
climb along the elemental arc at 89.5 degrees. 

Summarizing, the state incremental vector for an accelerating climb is 
given by 


m « 

AV 

= 

V. . ' v - 
l+l l 

Ah 


h. . - h. 
l+l i 

AY 


Y - Y 

i+1 i 

AR 


j[V i+ i cosYi + i + v i cosYi] 

AW 


j[W i+ l + Wi]At 

At 


RCi+x-Rq log f Rq ] 


6.4.1 Cruise Flight 

Cruise flight performance is computed by the Breguet equation. With constant 
velocity the distance travelled in time At is 



AR = V • At 

(151) 

Now 

SFC = | 

(152) 

so that 

_ AW 

~ (SFC)-T 

(153) 

Substituting in 

Equation (7.2.82) 



AR (SFC) -T AW 

(154) 

In cruise flight 

h~ w 

D T 

(155) 
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ad - v rh\ 

~ SFC V W 


On integrating 


! i,i - R i - spc = m lot <tt j 


C156) 


(157) 


Where the range factor W 7 is usually a slowly changing function of weight. 

NSEG.uses the inverse relationship to compute the weight change given 
a range increment 


r Ri+r R i, (R, rl -R.)-(SFC) 

L — RP J r— ili— i 1 

W. , * lie “ = Vhe *■ V(L/D) J 

l+l l 1 


(158) 


Alternatively, the program can be used with a time increment it by using the 
relationship 


r At* SFC, 

1 07 W 1 


W. . = W.e 
i+I i 

Several cruise modes are contained in the program including 

1. Constant altitude, constant Mach number cruise 

2. Constant altitude, constant C L cruise 

3. Constant Mach number, constant C^ 

Each of the three' cruise modes may be performed in the manner 

1. From Ri to Ri+i = ARi 

2. From to T i+1 = AT^ 

3. From Wj to Wi+i = AW^ 

A cruise flight is computed by summing over Ni steps. Thus, 

AR . = TaRi 

cruise £ 1 


(159) 


cruise 


■ ? 4T i 


cruise 


= l AW. 
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In all cases the total state increments are summed in the manner 


(AX) 


cruise 


lux). 


( 160 ) 


A mean range factor is used in all cruise calculations. The mean range 
factor, (RC^) , in each elemental arc bounded by {X}^ and {Xl^j is 
determined by an appropriate weighting of the range factors RCj_ and RC^+i 
which bound the arc. 


6.5 Descent 

The climb analysis of Section 6-4 is also used for the descent analysis. 
If the size of the flight path angle becomes too small (sinY < -1), the 
engine is throttled back to maintain a realistic flight path approximation. 


6.6 Level Flight Acceleration 

The approximate time to accelerate from Mi to Mi+l in level flight is 


t 


At. = a s (M i+1 - Mi)/Vi 
with a corresponding weight change 

(161) 

AW' = W. At' 
i 

so that 

' (162) 

W! , = W. - W. At' 
l+l l i 

Therefore, to the first order 

(163) 

v = g r T - D ] 

i+1 g o L W 1 

(164) 


The fuel flow at this point, Wj+i, can be obtained from the vehicle aero- 
dynamic and propulsion representation. 


V. 

l 


= — rv- 

7 1 1 


VP 


(165) 


and an improved estimate of the time to accelerate from Mj to Mi+j is 


At. = a (M. . - M. )/V. 
l s k 1+1 1 1 
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( 166 ) 



This gives an improved estimate of the weight change 


AH. = j{W i + 




At. 


and the corresponding range change 

4R i ■ K * W At i 

Summarizing, the level acceleration state increment is 


Ah 

= , 

0.0 

AV 


a sC M i+l " M i) 

AY 


0.0 



1 t * 

fil? 


jCKi * w i+1 ) it- 

5R 


* Mi) iti 
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<-> 
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»>iU - ♦ tT^ w l 
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(167) 


( 168 ) 


(169) 
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7. MISSION SEGMENTS 


The state incremental methods of Section 6 are used to create a variety 
of optional mission segments in NSEG. Each available mission segment 
option is briefly described below. All mission climbs, cruises, acceler- 
ations, and decelerations may be performed in forward or reverse direction. 
Each mission segment described below is performed as a distinct option in 
NSEG. There is some degree of overlapping capability in the available 
mission options. The mission option within NSEG is indicated for each 
mission segment for reference purposes. 


7.1 Generalized Climb or Descent (Mission Option 1) 

The generalized climb descent option incorporates several alternative 
flight path estimation methods. These methods include: 

Linear M-h path 

Constant dynamic pressure path 
Approximate optimal paths 

Each of the available paths is described below. 
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.7.1.1 Linear Climb 


This option climbs linearily from (M^, h.) to (M 2 , 1^) using a specified 
number- of linear cl^imb steps from (M^, hi) to (Mi + i, hi+i). The path is 
illustrated in Figure 14. 


7.1.2 Climb at Specified Dynamic Pressure 

This option climbs along a specified dynamic pressure line from (Mj hi) 
to (M 2 h 2 ] with appropriate terminal maneuvers. Along the constant 
dynamic pressure line a specified number of linear Mach-altitude segments 
are flown. Appropriate initial and final maneuvers are used when (Mi hi) 
or (M 2 h2) do not lie on the specified dynamic pressure line. The user 
may specify a climb at the terminal end point dynamic pressure. In this 
case, the final maneuver is not required. This path is illustrated in 
Figure 15. 


7.1.3 Rutowski Climb 

The Rutowski climb, Reference 1 flies from (Mi hi] to (M 2 I 12 ) along the 
path which must rapidly build up specific energy. If either of the points 
(Mi hi) and 'CM 2 h2) do not lie on this path, an appropriate terminal 
maneuver is employed. The Rutowski path is found by the following procedure. 

1. Compute the initial point specific energy 

Ei = Vi 2 /2g + h x (170) 

and find specific energy 

E 2 = V2 2 /2g + h 2 (171) 

and divide the energy change (E 2 - Ei) into N equal increments 

2. Search at each incremental energy level 

Ei = Ei + i • AE i = 1, 2, . . ., N (172) 

to find the point of maximum specific energy derivative, 

(Mi, hi) where 

^ = (Ti -Di) Vi/Hi (173) 

The E calculation is carried out for specified weight and 
load factor. 
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3. Fly a sequence of linear Mach-altitude flight increments 
joining the point (M^ hp and (Mi+i hi+j) 

A typical Rutowski path obtained from the program is illustrated in Figure 
16. The initial acquisition of fhe Rutowski path at 

h = h 1 + Ah ( 174 ) 

takes a vehicle from its initial condition to the Rutowski path with a 
velocity loss if this is required. The final maneuver may be either a 
transfer along a constant energy line from the Rutowski point at the final 
energy to the point M2I12. Alternatively, an altitude limit may be placed 
on the path such that when a Rutowski point lies above the final point, a 
transfer to the final point M2h2 occurs. These terminal maneuvers are 
sketched below. 





t> M 

(a) Initial Maneuvers 



The Rutowski path will observe both CL max and maximum dynamic pressure con- 
straints at the user's option. The thrust levels, vehicle weight, and load 
factors employed in the E calculation are specified by the user. Further 
details of this mission segment option may be found in Reference 2. 
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7.1.4 Maximum Rate of Climb 


A maximum rate of climb path between Mi hi and M 2 I 12 is generated in a 
similar manner to the Rutowski path of Section 7.1.3. However, in the 
maximum rate of climb path the search for maximum E is carried out at the 
constant altitudes 


h = + i * Ah i = 1, 2, . . ., N (175) 

where the altitude differential (h 2 -hi) has been divided into N equal 
increments. A typical maximum rate of climb path is shown in Figure 
16 . 


7.1.5 Maximum Acceleration 

A maximum acceleration path between Ml hi and M 2 h 2 is generated in a similar 
manner to the Rutowski path of Section 7.1.3. However, in the maximum 
acceleration path the search for maximum E is carried out at the constant 
altitudes 

M * Mj + i • AM i = 1, 2, . . N (176) 

where the Mach number differential (M 2 - Mi) has been divided into N equal 
increments. A typical maximum acceleration path is shown in Figure 
16. The maximum acceleration path satisfies the CL max and maximum dynamic 
pressure constraints of the Rutowski path. In addition, the condition 

AE i+ i * AE. (177) 

is imposed. That is, the sequence of points, Mj_ hi used in the acceleration 
will never produce a loss of specific energy. This is illustrated in 
Figure 17. 


'7. 1 . 6 Minimum Fuel Paths 

Minimum fuel path for given energy, altitude, and Mach number are obtained 
in a'manner similar to Sections 7.1.3 through 7.1.5, respectively. However, 
the search optimization criteria on E is replaced by the criteria 

<J> = Maximum [E/M] = = Max[^|] (178) 

When the search is carried out along lines of constant energy , the minimum 
fuel energy build up is found. YJhen the search occurs at constant altitude, 
the minimum fuel climb is found. When the search occurs at constant Mach 
number the minimum fuel acceleration is found. All appropriate terminal 
maneuvers and constraints described in Sections 7.1.3 to 7.1.5 are included 
in the minimum fuel paths . 
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7.1.7 Maximum Range Glide 


Hie maximum range glide path is obtained when the vehicle flies along 
the laws of the L/D contours tangency points to an appropriate path gen- 
erating surface such as constant energy, constant altitude, or constant 
Mach number. The maneuvers are thus similar to those of Sections 7.1.3 
to 7.1.5 using the optimization criteria 

4> = Maximum [L/D] (179] 

When the search is carried out along lines of constant energy, the maximum 
range glide for a given energy loss is found, Reference 1. When the search 
occurs at constant altitude, the maximum range glide for a given altitude 
loss is found. When the search occurs at constant Mach number, the maximum 
range glide for a given velocity loss is found. Some typical paths obtained 
from the NSEG program are illustrated in Figure IS. 


7.1.8 Range Biased Ascents 

Range biased ascents can be obtained when the vehicle flies along the locus 
of the T/(L-D) contours tangency points to an appropriate path generating 
surface. This can be seen as follows: 


and 


Now 


E = h + V 2 /2g 
m(~r) = T - D - W sinY 

-i-dE 


There from Equation (ISO) 

_ r V cosY • dE 
J dh V . dV 
dt + g dt 


cosY • dE 

1 dV 
sinY + — -rr 
g dt 


But from Equation 


1 

g 


(181) 


f dV, T-D 

l dt J W 


siny 


(180) 

(181) 

(182) 


(183) 


(184) 
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So that 


I 


R - / W C °| I ' 116 (185) 

Assuming that range biased ascents occur at small flight path angles with 
L - W, Equation (188) becomes 

R “ / ^ dE ( 186) 

Therefore, an energy-like approximation for a range biased ascent when 
(L/ (T-D) is a maximum at each energy level. It should be noted that when 
T - D = 0, no energy gain is possible; therefore, this singular condition 
must be avoided. In NSEG the per cent excess of thrust over drag which 
is acceptable is a program input. 

In a manner similar to Sections. 7.1.3 to 7.1.5 a range biased ascent 
between two energy levels occurs when the points of tangency between constant 
energy and T/(L-D) contours is flown. A range biased climb between two 
altitudes will fly the points of tangency between constant altitude and 
constant T/(L-D) contours. A range biased acceleration will fly the points 
of tangency between constant Mach number and constant T/(L-D) contours. 


7.1.9 Range Biased Ascents Based On 
Range Factor 

A second series of range biased ascents can be found on the basis of the 
range factor contours. These ascents are similar to those of Section 7.1.8 
with range factor replacing T/ (L-D) . 


7.2 Maximum Lift Coefficient 
Climb or Descent (Mission Option 2) 

The maximum lift coefficient path climbs from Mi hi to M 2 h 2 in N increments 
of altitude 

hi = hi + i • Ah i * 1, 2, ...» N (187) 

At each altitude the Mach number for maximum rate of climb using the angle 
of attack for ■* -S ^ oun< ^ 

M i - m *max rc (188) 

* 

The vehicle uses the linear Mach-altitude path path follower to climb between 
Mihi and Mi+i hi + i. 

Descents follow the same procedure as climbs, but in reverse order. 
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7.3 Radius Adjustment (Mission Option 3) 


This mission segment option performs an iteration on the range of one 
cruise segment to make the total range over the combined mission segments 
Sji j = Jl» J 2 • • • » equal to the total range over the combined 
mission segments S^; k = Kj, K 2 , . . . , Kjj. That is 


R = 



(189) 


where one of the ARj , and only one , is being adjusted to satisfy the range 
equality. 


7.4 Cruise Climb to Specified Weight (Mission Option 4) 

As an aircraft cruises at the Mach number and altitude for maximum range 
factor. Equation (151), the weight reduces. As the weight changes, the 
altitude for best range factor changes while the Mach number remains approx- 
imately constant. The altitude change results from the requirement to 
maintain the angle of attack for maximum lift/drag ratio. Thus, as the 
cruise progresses the altitude increases. 

The cruise may be performed in one step or it may be reduced to a sequence 
of five steps between latter case W = Wi and W = W 2 , Section 6.4.1. 

At the start of the i^h segment in this 

= W x + i • AW i = 1, 2, . . ., 5 (190) 

Each segment is flown at constant Mach number and lift coefficient and, hence, 
involves a climbing cruise. At the beginning of each cruise step the weight 
is instantaneously adjusted to the best altitude for the current weight. 


7.5 Cruise Climb for Specified Distance or Time (Mission Option 5) 

This mission segment option performs a cruise climb. Section 6.4 

for specified distance or time. The cruise may be performed with or without 

range credit. This form of cruise flight is performed in one step. 


7.6 Constant Altitude Cruise Between 
Two Weights, (Mission Option 6) 

This mission segment performs either 

1. Constant altitude, constant Mach number cruise 

2. Constant altitude, constant lift coefficient cruise 
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between two weights Wj and W 2 . The cruise is performed in one step, 
see Section 6.4 


7.7 Constant Altitude Cruise for 
Given Distance (Mission Option 7) 

This mission segment performs either 


1. Constant altitude, constant Mach number cruise 

2. Constant altitude, constant lift coefficient cruise 

between two distances Ri and R 2 . The cruise is performed in one step, 
see Section 6.4 


7 . 8 Constant Altitude Cruise for 
Given Time (Mission Option 8) 

This mission segment performs either 

1. Constant altitude, constant Mach number cruise 

2. Constant altitude, constant lift coefficient cruise 

between two times T. and T 2 . The cruise is performed in one step, see 
Section 6.4. This segment may be performed with or without range 

credit. 


7.9 Buddy Refuel Cruise (Mission Option 9) 

This mission segment determines the optimum in-flight refuelling point and 
how much fuel will be transferred. The tanker fuel off load capability is 
specified at three range/fuel combinations and a parabolic variation in 
available fuel as a function of range is assumed. That is, 

= a + bR + cR? (191) 

Cruise flight is assumed in any one of the three forms 

1. Constant Mach number, constant lift coefficient 

2. Constant Mach number, constant altitude cruise 

3. Constant lift coefficient, constant altitude cruise 

A maximum range for refuelling may be specified. Refuelling will occur at 
any point on the segment where 

1. Fuel receivable is greater than or equal to fuel available 

W Fr > Wp A (192) 
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2 . 


1 


Distance flown is equal to maximum refuelling range 

3. Minimum in-flight weight of the vehicle receiving fuel is 
reached where 

W MIN = W OWE + W PL + W p ' k F (193) 

where kp is the unusable residual fuel in the non-payload fuel. 
For refuelling purposes the maximum weight is taken to be the 
thke-off weight 


W MAX = W T0 


(194) 


7,20 Mach-Altitude-Weight Transfer (Mission Option 10) 

This mission segment option retrieves state components at the end of flight 
segment i and makes them available as the initial conditions for flight 
segment j. The initial conditions for segment j are thus a linear trans- 
formation of the final condition of segment i, 

{X}. = [P] i;j {X>. (195) 

Currently, the NSEG program is limited to a simple state component transfer 
on any combination of. the three components: Mach number, altitude, or 

weight. 


7.11 Alternate Mission Selection Option 
(Mission Option 11) 

This mission option retrieves either of two mission segments on the basis 
of terminal Mach number,' altitude or weight. Retrieval criteria may be 
based on any one of six possibilities: 


♦ = Min [Mj, M 2 ] (196) 

* 

<J> = Max [Mi, M 2 ] • ( 197 ) 

cf> = Min [hi, I 12 ] (198) 

4> = Max [hi, h 2 ] (199) 

4 > = Min [Wi, W 2 ] ( 20 °) 

<J> = Max [Wi, W 2 ] ( 2 °1) 


The segment to be retained is the one which satisfies the selected performance 
criteria. 
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7.12 Instantaneous Weight Change (Mission Option 12) 


& 


This mission segment option permits - an instantaneous change in vehicle 
weight, AW. The operation 

W i+ i - W. - AW (202) 

is performed. 


7.13 Instantaneous Mach/Altitude Change (Mission Option 13) 

This mission segment option provides an instantaneous change in vehicle Mach 
number, AM, and an instantaneous altitude change, Ah. The new Mach number, 
Mi+i, and altitude, hi+i, are specified directly; thus 


AM = Mi+i - Mj 
Ah = h^ + j - hj 


(203) 

(204) 


7,14 General Purpose and Point Condition Calculation 
(Mission Option 14) 

This mission segment option provides any of a variety of calculations 
described below: 

1. Best cruise altitude for given Mach number and weight 
based on range factor 

Max [RF ; M, W] (205) 

h 

2. Ceiling for a specified rate of climb at given Mach 
number and weight 

Max[RC; M, W] (206) 

h 

3. Mach number for maximum lift coefficient at given weight 
and altitude 

Max[C L ; W, h] (207) 

M 

4. Mach number for specified lift coefficient given weight 
and altitude 

Find [C^ ; W, h] (208) 

M 
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5. Maximum endurance Mach number given altitude, weight, and 
maximum lift coefficient 

Min[W; W, h, Cj^^] (209) 

6. Maximum Mach number at given weight and altitude 

Max[M; W, h] (210) 

M 

7 a. Mash number for maximum rate of clirrb at given weight and 
altitude 


Max[RC; W, h] 
M 


( 211 ) 


lb. Mach number for maximum rate of climb per pound of fuel at 
given weight and altitude 

Max[dh/dW; W, h] ’ r 


8. Approximate Mach number and altitude for maximum- range factor 
given weight 

Max[RF; W] ' (213) 

(M,h) 

9. Mach number for maximum range factor given altitude and weight 

MaxfRF; h, W] (214) 

M 

10. Various energy maneuverability parameters at specified load 
factor given Mach, altitude, and weight 

a. The required lift coefficient 

b. Specific excess power 

P g = E = (T-D)V/W 

c. Specific excess power divided by fuel flow 

P S /W = E/W = (T-D)V/ (W • W) (215) 

d. Specific excess power divided by fuel flow and multiplied 
by fuel remaining (AE capability) measure 


(216) 


W F E W F 

p s — = -r~ 
w w 


dE 

dW 


W p « AE 


e. Specific energy 

E s = h + V 2 /2g 

f. Load factor at Pg=0.0 

g. Steady state turn radius computed as follows: 

= C^, for given load factor 

Now for given bank angle, B A 

W = qSC^ • cos(B a ) 

and the centrifugal force is 

Lsin(B A ) = WV^ _ Lcos(B A )V 2 

Rg ~ Rg 

V 2 

*** R ~ gtan-B A 
but from Equation (7.2.149) 


cosB. = - 

A gSC L 


tanB A =“^(5|Et) 


1.0 


Substituting Equation (7.2.153) and (7.2.151) 


(217) 


(218) 


(219) 

( 220 ) 


( 221 ) 


( 222 ) 


V 2 / 1.0 

R - r\ (223) 

It should be noted that this mission segment option may employ directly 
specified values of Mach number, altitude, and weight of these state components 

may be picked up from the previous mission segment termination. The option 
to reset Mach number, altitude, and weight from any previous segment termi- 
nation is also available within the option. 
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7.15 Iteration to Fly a Specified Distance (Mission Option 15) 

This mission segment option perturbs the range increment in segment i to 
provide a specified total range (from mission initiation) in segment j 



This is illustrated above where AR^ is perturbed to satisfy the condition 


within an error of one nautical mile. 


7,16 Climb or Accelerate (Mission Option 16) 

This mission segment option provides a climb or acceleration between two 
Mach number-altitude points (Mi hi) and (M2 ^2) • These two flight conditions 
must be defined in two mission segments, segment i and segment j. The 
climb or acceleration will then join the two points. Climb or acceleration 
paths may be performed in either a forward or reverse time direction. Descents 
are not permitted. The mission segment option may be performed with or without 
range credit. 

Fuel burning decisions are made according to Mil-C rules while going from 
condition 1 to 2. Thus, fuel is burned if 
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7.17 Fuel Weight Change (Mission Option 17) 

A computed or specified fuel weight change is introduced through this 
mission segment option. The operation performed is 



W i+ 1 = V. - AW 

(224) 

T . +1 = T. + AT 

The option can be used to compute 

(225) 

1 . 

Loiter fuel requirements 


2. 

Warm up and take-off fuel 


3. 

Combat fuel 



Take-off fuel when computed is carried out through the take-off analysis of Section 
6.1. If a detailed take-off analysis is not required the option of Section 

7.18 is used. Warm up fuel calculation is computed for given power 
setting and time. Loiter* fuel calculation is for flight at specified Mach 
number, altitude, weight, and a given time. Combat fuel calculation is for 
specified time or degrees of turn at a given load factor. If the degree of 
turn option is used, the following calculation is performed. 

L = n • W (226) 

where n is the load factor. The centrifugal force is 

f r =Vl 2 - w 2 (227) 

and the turn radius is 

R = WV 2 

g/L2-W2 (228) 

The thrust force is set to drag at the turn 

T = D (229) 


7.18 Fuel Allowance (Mission Option 20) 

This mission segment option computes the fuel allowance for a specified time 
at 

1. Given power setting 

2. Given thrust/weight 


7.19 Engine Scaling (Mission Option 21) 

When using the turbojet, ramjet and/or combined engine option of Section 2.2, 
the engines may be resized at any point in the mission by specifying this 
option. 


7.20 Generalized Iteration Control Option (Mission Option 22) 

Option 22 can be used to provide an automatic iteration loop within the total 
mission starting and finishing with arbitrary segments input in Option 22 
segment data. The purpose of the iteration is to vary a specified quantity 
in the starting segment until a specified quantity in the final segment 
specified achieves a desired value. 



7.21- Thrust Specification in Mission Segment Options 




The vehicle propulsive representations have been discussed in Section 2. 

In the first option of Section 2.1 there are three available maximum thrust 

tables, T . ; j = 1, 2, 3. These tables are referenced as follows: 
maxj J 


T , = maximum dry thrust 
maxi J 

T _ = maximum wet thrust 
max2 

\iax3 * max ^ mum P ower 


C230) 


Throttling may only be used for T max i and T max2 - In using the various mission 
segment options an appropriate choice of thrust must be made. The options are 


1. T = D 

2. T = maximum dry 

3. T = maximum wet 

4. T = maximum power 

5. T = thrust for given power setting, dry. 

In the second option of Section 2,2 either throttled or maximum 
used. 


(231) 


power may be 


8.0 FLIGHT ENVELOPE CALCULATIONS 

Several gross flight envelope calculations may be performed. All flight 
envelope computations are subject to the conditions 

C L ^ ^ L lim ' coefficient limit 

M < Miim , Mach number limit 
q < q lim , dynamic pressure limit. 

Propulsive and aerodynamic characteristics must be specified. 


8.1 Climb Path History 

Given an initial weight, warm up, and take-off fuel allowance, a maximum 
rate of climb path is performed from 

P 1 = p iQ*i» 0.0) (232) 
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to 


P 2 = P 2 CM 2‘ RF^ C 233 ) 

where 

^MAX RF = alt ^ tu< ^ e f° r best range factor at M 2 
Alternatively, P 2 may be selected as 

P 2 = ^“mAX RF» h MAX RF^ C 234 ) 

The calculation is performed in ten equal altitude increments from Pj to 
P 2 - Climb paths are generated for N distinct weights 

Wi = W 0 + i • AW; i - 0, 1, .... N-l ( 235) 

8.2 Endurance versus Weight at Various Altitudes 

The endurance is calculated at a given altitude for the weights W^ = W 0 +i*AW; 

i f 1, 1, . . . , N-l. • Mach number selected is for best endurance. 

The calculation may be repeated for any number of altitudes; h = h 0 +i*Ah; 
i = 0 , 1, . . . . 

8.3 Optimum Cruise Climb at Various Mach Numbers 

An optimum cruise climb between Wj and W 2 in a specified number of weight 
increments. The path is repeated for an array of Mach numbers and altitudes 

M i ■ M o ♦ i • iM - ^ 0 , 1 , 2 , C256) 

bj = h Q + j • Ah; j = 0, 1, 2, (237) 
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8.4 Contour Presentation Capabilities 


A set of point calculations (vehicle capability at given flight conditions) 
are carried out over a two-dimensional array of Mach-altitudes, Mj, hj . 

The resulting matrix of capabilities,, F^j , is then supplied automatically to the 
NSEG contour plotting routine CONPLOT, and the contours of the function F k 
in the Mach-altitude plane are obtained in the form of CALCOMP, Houston 
plotter, or CRT display device output. At the present time twelve functions, 

Fl to Fl2, may be output in contour form. Each contour plot is described 
briefly below. 

8.4.1 Specific Energy Time Derivative, E, (INDMAP=1) 

The specific energy time derivative is computed according to the expression 

F 1 = E(M, h) - (T - D)V/W (238) 


where 

E = energy total time derivative 

T * thrust obtained at a specified power setting or at T * D; 
wet, dry, or maximum power options are available 

D = drag computed for a specified load factor 

V = flight velocity 

W = aircraft weight 

Some typical energy derivative contours for a large four-engine transport are 
presented in Figure 19. The minimum contour shown is for the condition 
T - D = 0. Hence, the flight envelope is a by-product of the E map when 
suitable constraints such as , and dynamic pressure limits are added. 


8.4.2 Specific Energy/Fuel Flow, E/m, (INDMAP=2) 


The E/m contour presents the specific energy time derivative over the fuel 
flow map. Since 


A/1 _ dE/dt _ dE 
E/m = dS7dt = dm 


(239) 


The map illustrates an aircraft's ability to convert fuel into energy at 
specified flight conditions. 


The point calculation employed is 

F 2 = E/m = (T - D)V/(Wm) 


(240) 
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where m is the fuel flow rate. The various thrust and drag options discussed 
in Sections 2 and 3 may be employed to produce a family of maps. A 
typical example for the large subsonic transport at maximum thrust and lg 
flight is shown in Figure 20. 

8 . 4.3 Lift/Drag, L/D, (INDMAP=3) 

Lift/drag contours present a measure of the airplane's aerodynamic efficiency. 
The L/D maps indicate its range capability in unpowered flight and partially 
reflect the cruise range capability. Mass can be produced for any specified 
load factor. A typical contour for the large subsonic transport in level 
flight is presented in Figure 21. 

F 3 = L/D (241) 


8.4.4 Range Factor, Rp, (INDMAP=4) 

Range factor contours present a measure of vehicle cruise range capability. 
Maps are produced for level flight with thrust equal to drag at a specified 
aircraft weight. 

F 4 - Rp - C^c) ( 5 ) (242) 

where SFC is the specific fuel consumption. The user may elect to construct 
maps for other than level unaccelerated flight. However, the interpretation 
of these charts is not clear. A typical unaccelerated flight range factor 
contour map for the large subsonic aircraft is presented in Figure 22. 

8 . 4.5 Thrust (INDMAP=5) 


The thrust map is available as a device for examining thrust input data or 
the thrust component of other mapped functions . The map can be obtained for 
wet, dry, maximum, or throttled power setting. The maximum power thrust 
map for the large subsonic transport is presented in Figure 23 . 


F 5 = 


wet 1 dry’ 


max 


or T, 


PS 


(243) 


8.4.6 Drag Map (INDMAP= 6 ) 

The drag map provides a device for inspecting, drag data, input or the drag 
component of any other map. Drag maps are produced for a specified load 
factor. A lg drag map for the large subsonic transport is presented in 
Figure 24. 

F 6 = D g = 1 (244) 

8.4.7 Specific Fuel Consumption, SFC, (INDMAP=7) 


Specific fuel consumption maps are provided as a data input inspection device 
or as an aid to visualizing the specific fuel consumption component of other 



maps. Maps may be obtained for wet, dry, maximum, or throttled power settings. 
Maximum power specific fuel consumption of the large subsonic transport is 
presented in Figure 25. 

F 7 = SFC wet ’ SFC dry’ SFC max’ ° r SFC PS (245) 

8.4.8Fuel Flow Rate, m, (INDMAP=8) 

The fuel flow maps are provided as a data input inspection device or as an 
aid to visualizing the fuel flow component of other maps. Maps may be obtained 
for wet, dry, maximum or throttled power settings. Maximum power fuel flow 
for the large subsonic transport in level unaccelerated' flight is presented 
in Figure 26-. 


8.4.9 Specific Energy (INDMAP=9) 

The specific energy map 

F 9 = E = h + V 2 /2g (246) 

is provided as a user's convenience in visualizing the trajectory points 
between constant energy lines and any other set of contours. An example is 
presented in Figure 27. 


8.4.10 Lift/ (Thrust - Drag), L/(T-D) (INDMAP=10) 

The lift/ (thrust - drag) contours are useful for determination of maximum range 
powered flight. 

Assuming that maximum range flight occurs at small flight path angles 

R “ / tCd dE (247) 

Therefore, the energy- like approximation to maximum range flight occurs when 
L/ (T-D) is a maximum at each energy level. It should be noted that when 
T - D - 0, no energy gain is possible; therefore, this singular condition must 
be avoided. In NSEG the per cent excess of thrust over drag which is 
acceptable is a program input. A typical L/(T-D) map for the large subsonic 
transport is presented in Figure 28. 

F - (L/CT-D) S p ecified excess thrust (248) 

8.4.11 Turn Radius (INDMAP-11) 

Turn radius maps give a gross indication of aircraft's combat capability. 

Turn radius is computed by equating the aircraft's lift capability in steady 
state of decelerating flight using the following expression 



(249) 



where Cl is determined so that (a) thrust equals drag for steady state 
flight and (b) Cl equals Cl maximum for minimum instantaneous turn radius. 

Typical radius of turn map for the subsonic transport are presented in Figure 

29. 


8.4.12 Time to Turn (INDMAP=12) 

Time to turn through 180 degrees is presented as a supplement to the turn 
radius map. When the minimum instantaneous turn radius calculation is employed, 
the maps do not give a true time to turn. They merely indicate how long a 
time the aircraft would take to turn if it could maintain its current turn rate. 
When steady state turns are considered, true time to turn is obtained which 
will frequently be much longer than is required for a decelerating turn. 

Typical time to turn maps for the subsonic transport are illustrated in Figure 30. 
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Hi-Lo-Lo-Hi Radius Mission 



1 . Take-off 

2. Maximum rate of climb to best cruise altitude given weight and Mach number 

3. Constant Cl climb to best cruise altitude for new weight and Mach number 

4. Breguet cruise to given range, R 

5. Instantaneous state change to dash Mach number and altitude 

6. Constant Mach number-altitude cruise to total range, R 

7. Drop ordnance, instantaneous weight change 

8. Constant Mach number- altitude return cruise to given weight 

9. Maximum rate of climb to given Mach number-altitude 

10. Breguet cruise to given weight 

11. Instantaneous state change to best endurance Mach number for given altitude 
and weight 

12. Loiter for given time 


FIGURE 1. TYPICAL NSEC 


MISSION PROFILE 



ALTITUDE 



FIGURE 2. TYPICAL FLIGHT POINT PERFORMANCE MAPS 
(SPECIFIC ENERGY AT VARIOUS LOAD FACTORS, WEIGHTS,. AND THRUSTS) 
WITH TRAJECTORY LIMITS OVERLAID 
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FREESTREAM DIRECTION 



FIGURE 3. REFERENCE STATIONS AND GEOMETRY 






PARALLEL TO 
WING UNDERSURFACE 


= SPILLED AIRFLOW/P 1 V 1 


2 = Aj^/Tan Sj. 

1 = ( P 2 - V A 1 

'SPILL = F 2 COSy - F l Siny 


FIGURE 4. SPILLAGE LIFT 
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PARALLEL TO 
WING UNDERSURFACE, 
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PILL = F 1 COSV + F 2 SinV 


DRAG COMPUTATIONS 





















Turbojet Fuel Flowrate 
Air Flowrate Ratio 


TABLE VII 


TABLE VIII 
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\A " K07 A * CO (3-DIM.) 


A - .58 


A*0 (THEORY) 


FIGURE 9. TYPICAL FIRST PEAKS IN LIFT COEFFICIENT VERSUS ANGLE 

OF ATTACK 
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FIGURE 10. SUBSONIC MAXIMUM LIFT OF WINGS WITH POSITION OF 
MAXIMUM THICKNESS BETWEEN 35 AND 50 PER CENT CHORD 
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Figure 14. Linear Mach Altitude Mach Segment 
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Figure 17. Energy Constraint Imposed on Maximum Acceleration Path 
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W=575,000 lbs. 
LEVEL FLIGHT 





MACH NUMBER 

FIGURE 27. SPECIFIC ENERGY MAP (COMPUTER DRAWN) 



MACH NUMBER 


FIGURE 28. LIFT/ (THRUST-DRAG) (COMPUTER DRAWN) 



ALTITUDE x 10 , FT. ALTITUDE x 10 , FT. 
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